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Within a microscopic two band model of planes and chains with a pairing potential in the planes and 
off diagonal pairing between planes and chains we find that the chains make the largest contribution 
to the Josephson tunnelling current and that through them the d-wave part of the gap contributes to 
the current. This is contrary to the usual assumption that for a d-wave tetragonal superconductor 
the c-axis Josephson current for incoherent tunnelling into an s-wave superconductor is zero while 
that of a d-wave orthorhombic superconductor with a small s-wave component to its gap it is 
small but non-zero. Nevertheless it has been argued that the effect of twins in YBCO would lead 
to cancellation between pairs of twins and so the observation of a current in c-axis YBCO-Pb 
experiments is evidence against a d-wave type order parameter. We argue that both theory and 
experiment give evidence that the two twin orientations are not necessarily equally abundant and 
that the ratio of tunnelling currents in twinned and untwinned materials should be related to the 
relative abundance of the two twin orientations, 
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I. INTRODUCTION 

The symmetry of the order parameter in the cuprates 
is the subject of much current debate. While several 
experimental techniques such as ARPES and cits 

1^,^ probe the magnitude of the superconducting order 
parameter Josephson tunnelling is, perhaps, the only ob- 
servable phenomena which directly probes the phase of 
the superconducting order parameter as well as, indi- 
rectly, its magnitude. 

The focus of the debate has been upon whether the 
order parameter has s-wave (even under a 90° rotation) 
or d-wave (odd under a 90° rotation) symmetry. Since 
most experiments are performed upon orthorhombic ma- 
terials in which the s- and d-wave symmetries belong to 
the same irreducible representation there is no clear dis- 
tinction between them and, rather than the symmetry 
of the order parameter, the discussion should focus on 
whether there are nodes present in the order parameter 
which cross the Fermi surface. 

Although the s- and d-wave representations of the 
tetragonal systems mix freely in orthorhombic systems 
we can still speak of the "s" and "d" parts of the order 
parameter if what we mean are, respectively, the parts 
which are odd and even under a 90° rotation. It would, 
however, be highly unusual for one of these components 
to be present in an orthorhombic superconductor with- 
out the other being present as well. Further, if the odd, 
or d, part is dominant then nodes which cross the Fermi 
surface will be present in the order parameter (what we 
call a (d -f es)-wave order parameter) while if the even, 
or s, part is dominant the nodes, if they are present, will 
not cross the Fermi surface (what we call an (s-f ed)-wave 
order parameter). It is possible for an order parameter 
that is even under a 90° rotation (ie, s-wave) to have 
nodes but we consider this situation unlikely and we will 



not examine this possibility here. 

Several corner junction YBCO-Pb tunnelling experi- 
ments are cited as strong evidence for the order pa- 
rameter in YBCO having d-wave type symmetry. Other 
tunnelling experiments which contain only YBCO in a bl- 
and tri-crystal rings |l^-|l^ also give strong indications 
of order parameters with d-wave type symmetry. Both 
of these types of experiments have the plane of the junc- 
tions perpendicular to the a-b plane - what we refer to as 
"edge junction" experiments. It is the observation of tt 
phase shifts in the corner junctions and half flux quanta 
in the ring experiments that is the strong evidence for 
the order parameter in YBCO having d-wave type sym- 
metry, jl^ Note that for the observation of tt phase shifts 
there is no difference between a d-wave order parameter 
and a (d -|- es)-wave order parameter. 

The results of these edge junction experiments are in- 
dependent of the presence of twins in these materials, jlj] 
This indicates that the order parameter is phase-locked 
across the twin boundaries, although calculations indi- 
cate [|l^,|l^ that the magnitudes (as well as the relative 
magnitudes and perhaps the relative phase as well) of 
the d and s may change in the twin boundary. Although 
this phase-locking is not unexpected it may have other 
consequences. 

There is a second class of YBCO-Pb tunnelling ex- 
periments in which the plane of the junction is parallel 
to the a-b plane - what we refer to as "c-axis junction" 
experiments. The presence of a current in these 

experiments is cited as evidence against a d-wave order 
parameter in the literature | |l6| , p^ a s well as articles for 
more general physics audiences. |20[ In fact, due to the 
orthorhombic symmetry of these materials the order pa- 
rameter will have (d -I- es)-wave symmetry and a c-axis 
tunnelling current of reduced magnitude is expected in 
untwinned materials. ||2l|,p2| The current is not just due 
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to the es part of the order parameter but is also caused by 
the d part of the order parameter due to the orthorhom- 
bicity of the Fermi surface, a fact that appears not to be 
widely appreciated. 

While there is quite a large variation in the maximum 
resistance-tunnelling current product, IcRn, observed in 
these c-axis tunnelling experiments, IqRn in the un- 
twinned samples is, on average, about twice that in the 
twinned samples. |24[ The argument that a finite c-axis 
YBCO-Pb tunnelling current is evidence against a d-wave 
order parameter is based upon the intuitive assumption 
that in the presence of twins any orthorhombicity in a 
phase-locked order parameter will average out over the 
junction area. In fact, theory clearly shows that this 
would be an unusual condition and that one twin should 
be significantly more abundant that the other. An ex- 
periment which measures the relative abundance of the 
two twins by comparin g th e intensities of the (200) and 
(020) diffraction peaks ||2^ indicates that one of the two 
twins is, depending upon the specific sample, 2-3 times 
more abundant that the other. |^ 

If the ratio of the two possible twin orientations is mm 
and if the current-resistance product, IcRn{T — 0), is 

/c/J^Y^ for an untwinned junction then for the twinned 
sample the current should be: 



is the the pair susceptibility, with: 



IcR 



(twin) 
N 



IcR 



(1) 



This is consistent with the experimental observation 
of Sun et al. [p^ in which the Josephson tunnelling cur- 
rent in untwinned c-axis YBCO-Pb is about double that 
observed in twinned junctions if the ratio of twins is ap- 
proximately 3:1. 

In Section II we present a BCS formalism for a model 
system in which a tetragonal Cu02 plane layer is coupled 
to an orthorhombic CuO chain layer and give the equa- 
tion for the resistance-tunnelling current product, IqRn- 
In section III we present the results of some representa- 
tive calculations as well as some experimental results on 
the relative abundances of the two twin orientations. In 
section IV we make some concluding statements. 



II. FORMALISM 

For a bilayer system (ie, a = 1,2) the coupled BCS 
equations are: [g2[ 

^k,l = ^ X! (^k,q,llXq,l + Vk,q,12Xq,2) 



q 

^k,2 = ^ X/ (^k,q,12Xq,l + ^k,q,22Xq,2) 

q 

where we have taken Vk,q,i2 — Vk,q,2i and: 

Xq.ct = (OqT.a^— qi.ct) 

■ tanh 



(2) 
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where the band energies in the normal state. 

Note that we have taken (aq|_ia_q|^2) = 0, ie there is 
no pairing between electrons in different planes although 
other workers do not make this assumption. 

We note that this set of equations (^) is invariant under 
substitution {/ik,2, Vk,q,i2} — > {-^k,2, -14;,q,i2} which 
means that the overall sign of Vk,q,i2 only affects the rela- 
tive sign of the order parameters in the two layers and not 
their magnitudes. This is interesting because it means 
that the effect on Tc of having an interlayer interaction is 
independent of whether this interaction is attractive or 
repulsive, although the c-axis Josephson tunneling cur- 
rent still depends upon the relative sign of the interlayer 
interaction. 

The c-axis Josephson junction resistance-tunneling 
current product, IcRn{T), through a superconductor- 
insulator-superconductor junction for incoherent c-axis 
tunneling is given by the relation: p7|] 



^cnNKl) ^(Pb)(o)^(YBCO) (0)7^2 

x^A(P>^)(a;„)A(™)(c.„), 

TJ. 



where: 
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in which the superscript (•) indicates on which side of the 
junction the dispersion and order parameter are on, the 
sum over ujn = ■nT(2n — 1) is for all Matsubara frequen- 
cies, R is the resistance of the junction and N^'\0) is the 
normal state electronic DOS given by: 



lim 



r 



r'^'o ttQ ^ _ ^)2 + j^2 



Since the DOS and ^^^^^^ for lead are constant the sum 
in Eq. Hcan be performed and A^^^^ (a;„) is given by: |2^] 



(5) 



If the tunneling were coherent the matrix element 
(which is incorporated into i?) would have a (k — k') de- 
pendence, and the sums over k-space wouldn't be sepa- 
rable. 
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III. RESULTS 

We use the same band structure as in our previous 
work on the penetration depth of a coupled chain-plane 
bilayer in which we assume tetragonal symmetry for 
the electron dispersion in the Cu02 plane with first and 
second neighbour hopping, while the chains are quasi- 
one dimensional with very different hopping probabili- 
ties in the x- and y- directions (the chains are along the 
y-direction). For the pairing potential, Vk.q.Q/3, which 
appears in the coupled BCS equations (||) and which de- 
termines the superconductivity, we use a nearly antiferro- 
magnetic Fermi liquid model with magnetic susceptibil- 
ity of the phenomenological form given by Millis, Monien 
and Pines (mmp): ||2^ 

t/k,q,a/3-5a/3^^^2|k_q_Q|2' 

where is the magnetic coherence length and Q is the 
commensurate vector (tt, tt) in the 2-D Brillouin zone. 
Our results are not specific to this interaction - any in- 
teraction which would result in a d-wave order parameter 
in the Cu02 plane layer will give a, (d + es)-wave order 
parameter when the Cu02 planes are coupled to the CuO 
chain layers. 

No pairing interaction is assumed to act directly in 
the chain band, ie 322 = 0, so that the superconductiv- 
ity in the chains is entirely due to the gi2 — 1721 7^ 0. 
This parameter ((712) is fixed to get a critical temper- 
ature value of fOOK for a chosen value of the in plane 
pairing, gn. In Fig. |^, we show results for the gap 
value Z\k.i in the planes (Fig. |l|a) and for Z\k.2 in the 
chains (Fig. [^b) as a function of k in the first Brillouin 
zone. The order parameters result from a numerical so- 
lution of equations (^). For the runs shown in Fig. |^, 
{ffii, 312,522} = {26.2,10,0} and T, = lOOK. 

On the right hand side of Fig.|l|, we have decomposed 
the order parameters into d-wave (c) and (d) and s-wave 
(e) and (f) components for the plane and chain, respec- 
tively. While the main component is certainly the d-wave 
part, the s-wave admixture is, nevertheless, significant in 
magnitude in both the chain and plane bands. 

A useful representation of these gap results is to show 
the contours of gap zeros on the same plot as the Fermi 
surface. This is presented in the series of frames shown in 
Fig. 1^. The top frames apply to the plane while the bot- 
tom frames apply to the chains. In all cases, (a), (c), (e) 
for the planes and (b), (d), (f) for the chains, the same 
Fermi surface (dashed curves) was used. By choice, the 
Fermi contour have tetragonal symmetry in the top fig- 
ure while the chain Fermi surface is a quasi straight line 
along kx as is expected for chains along the y-direction in 
configuration space. The pictures are for three different 
values of pairing potential. The first set of two frames 
(a) and (b) are for {511,312,322} = {29.9,5,0}, i.e. very 
little coupling between chains and planes (off diagonal 
512 small). In this case, the gap in the plane is nearly 



pure d-wave as is also the induced gap in the chains. The 
zero gap contours are given by the solid line and would 
be the main diagonals in a pure d-wave system. As the 
coupling gi2 is increased {311,512,522} = {26.2, 10,0}, a 
significant s-wave component gets mixed into both solu- 
tions and the gap nodes move off the main diagonals of 
the Brillouin zone (This is the solution that is plotted in 
Fig. |l|). The gap nodes still cross the Fermi surfaces in 
both chains and planes. As the coupling is further in- 
creased to {511,312,522} = {9.18,20,0}, Fig. I (e) and 
(d), the gap nodes move far off the diagonal and for the 
chains they no longer cross the Fermi surface so that 
there is a finite minimum value of the gap on this sheet 
of the Fermi surface. 

Calculation of the product of Josephson current times 
resistivity given in equation (|) for a YBCO-I-Pb tunnel 
junction gives for the intermediate case {511,312,522} = 
{26.2, 10, 0} the tunnelling current, IcRn{T = 0), is 0.2- 
0.3 meV for the plane layers and ±(2-3) meV for the chain 
layers. The relative sign depends upon whether the in- 
terlayer interaction, 1712, is attractive or repulsive. These 
values for chain and plane tunnelling current, IcRn, are 
upper and lower bounds and note that the expected value 
should be somewhere in between since the IcRn{T = 0) 
values are dependent upon the amount of tunnel junc- 
tion area which is covering the uppermost chain or plane 
layer. 

The order of magnitude difference obtained between 
plane and chain layer for the IcRn{T = 0) product is 
understood as follows. For a given Matsubara frequency, 
LUn in equation (|^), the contribution to the sum over k, 
which ranges over the entire Brillouin zone, is strongly 
peaked about the Fermi surface because the denominator 
in (P becomes smallest in this case. This is seen clearly 
in Fig. H where we show the integrand of equation (|^), 
again for the intermedicate C&>S6, clS 3j function of k in the 
first Brillouin zone for planes (a) and chains (b), respec- 
tively for a particular Matsubara frequency, w„ =50meV. 
It is clear that in the plane layer the positive and negative 
parts largely cancel each other. They would, in fact, give 
exactly zero if the gap had pure d-wave symmetry. The 
situation for the chains is completely different because 
the Fermi surface now does not have tetragonal symme- 
try and even if the gap was pure d-wave there would not 
be a large degree of cancellation between the positive and 
negative regions. Thus, we note the important result that 
the main part of the Josephson current coming from the 
chains is due to the d-wave part of the gap function and 
would still be large if we did not account for the s-wave 
admixture. Thus, in an experiment on an untwinned sin- 
gle crystal of YBCO the Josephson current coming from 
the chain part of the Fermi surface is sampling mainly the 
d-wave part of the gap and therefore, such experiments do 
not reflect directly the s-wave admixture. 

It has been argued that for twinned samples, the c-axis 
Josephson current will cancel because of the cancellation 
between a- and 6-twins. This argument would apply 
equally well to our work, since we should then average 
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over pictures as shown in Fig.^(b) with opposite phases 
(a- and 6-twins). However, we do not expect twins to be 
present in equal numbers, and the expected cancellation 
will not occur, psj] 

High resolution x-ray scattering measurements have 
been carried out on a small single crystal of 
YBa2Cu3 6.93 for the purpose of investigating its twin 
structure. Measurements were made with an 18kW rotat- 
ing anode x-ray source and a double-axis diffractometer 
using a perfect Ge(lll) monochromator with sufficient 
resolution to cleanly separate Cu K^i from Cu Kq,2 radi- 
ation. Full details will be given in a separate publication. 

The crystal under study was grown by the UBC group 
and had approximate dimensions of 1 x 1 x 0.03 mm"^ , with 
the c-axis oriented along the thin dimension. The crys- 
tal was sufficiently thin that scattering within the or- 
thorhombic basal plane (such as (2,0,0) and (0,2,0)) could 
be performed in transmission geometry and such mea- 
surements thus probe the bulk of the crystal. Represen- 
tative longitudinal scans of the (2,0,0) and (0,2,0) Bragg 
peaks are shown for both twin domains in Figure 4. As is 
clear from this figure, one twin domain has a peak inten- 
sity which is more than an order of magnitude stronger 
than the other. This leads naturally to an interpreta- 
tion of the twin structure in the crystal in terms of a 
majority and a minority domain. Also, as would be ex- 
pected in this scenario, the minority phase lineshape is 
very noticeably broader than that of the majority phase, 
whose lineshape appears to be approximately resolution 
limited. This is shown in the inset to Figure |^, where the 
same data is plotted on a linear intensity axis, and has 
been scaled so that the peak intensity of the two domains 
coincide. 

The majority and minority phase domain distribution 
in this crystal was found to be inhomogeneous. X-ray 
scattering measurements were performed with a very nar- 
row (~0.05 mm) incident beam to allow measurements 
which probe different regions of the crystal. The twin 
structure was investigated as the sample was translated 
through the narrow beam along the majority phase a- 
axis. This investigation yielded results which ranged 
from completely untwinned on one extreme edge of the 
crystal to an approximate 1 to 1 ratio of the volume frac- 
tion of majority to minority domains on the other edge. 
Interior regions of the crystal yielded some intermedi- 
ate value of this ratio. We estimate the average ratio 
of the volume fraction occupied by the majority and mi- 
nority twin domains for the entire sample to be about 
2-3 to 1. While these measurements reveal a complex 
inhomogeneous morphology to the twin structure, it is 
certainly clear that macroscopic regions of the crystal 
exist in which one domain predominates over the other. 



IV. CONCLUSIONS 

For c-axis incoherent YBCO-TPb Josephson tunnelling 
junctions, the d-wave component of the gap parame- 
ter can contribute very significantly because of the or- 
thorhombic nature of the chain Fermi surface which em- 
phasizes the contributions to the Josephson current from 
those parts of the Brillouin zone along the Fermi surface. 
The s-component of the gap in the chains and planes will 
also contribute but this may be less important so that 
such experiments on untwinned samples do not probe di- 
rectly the s-wave admixture of a predominantly d-wave 
gap function. 

We argue that the observation of a c-axis tunnelling 
current in the experiments of Sun et al. p?! and Tanaka et 
al. p8[ | on twinned samples can also be understood within 
our theory because we do not expect equal numbers of a- 
and 6-twins in their experiments. 
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FIG. 2. The Fermi surface (dashed curves) and gap nodes 
(solid curves) for a Cu02 plane layer (top frames) and a CuO 
chain layer (bottom frames) for three different interlayer in- 
teraction strengths (left, middle and right frames). As the 
strength of the interlayer interaction, gi2, is increased (left to 
right frames) the proportion of the s-component in both layers 
increases. If the interlayer interaction were further increased 
the gap nodes would leave the Brillouin zone altogether and 
the order parameter would become s-like. Note that the Fermi 
surface in the Cu02 plane layer is tetragonal but that the gap 
node is orthorhombic. 




FIG. 1. The order parameter in the first Brillouin zone 
for (a) the plane layers and (b) the chain layers. Beside are 
the projections of the d components, (c) and (d), and the s 
components, (e) and (f). The vertical scale in all frames is the 
same. Note that the relative phase of the d-components in the 
two layers are the same while that of the s-components are 
opposite; this is caused by the interlayer interaction, Vk,q,i2, 
being negative (ie, repulsive). 




FIG. 3. The k-space integrand of ^(ybco) ^^^-^ ^^^g |^ 
for (a) the Cu02 plane layer and (b) the CuO chain layer. In 
the plane layer the positive and negative parts mostly cancel 
and the resulting current is small while in the chain layer 
only a very small amount of the integrand, A'^^^^'' (tj„), is 
negative and the resulting current is large. This effect is due 
to the different Fermi surfaces in the two layers. 
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FIG. 4. Representative longitudinal x-ray scattering 
scans through the (2,0,0) and (0,2,0) Bragg peaks of the 
YBa2Cu3 06.93 single crystal are shown on a logarithmic scale. 
It is clear that the peak intensity of the majority phase twin 
domains arc at least an order of magnitude stronger than that 
of the corresponding minority phase domain. The inset shows 
the same data, but now on a linear intensity axis, and scaled 
so that the peak intensities for the two domains coincide. The 
Bragg features of the minority twin domains are clearly much 
broader than those of the majority phase domains. 
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